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Abstract. For a functional £ and a peak selection that picks up a global maximum of 
S on varying cones, we study the convergence up to a subsequence to a critical point 
of the sequence generated by a mountain pass type algorithm. Moreover, by carefully 
choosing stepsizes, we establish the convergence of the whole sequence under a "local- 
ization" assumption on the critical point. We illustrate our results with two problems: 
an indefinite Schrodinger equation and a superlinear Schrodinger system. 
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1. Introduction 

Let us consider 3V a Hilbert space with inner product ( • | ■ ) and norm 1 1 • 1 1 , and a functional 
& £ If (Jf;M.). In this work, we develop a provably convergent "general" mountain pass 
type algorithm to approximate saddle points of S ', with a Morse index possibly larger than 
one. The pioneer work in this direction is due to Y. S. Choi and P. J. McKenna [4] who 
proposed a constrained steepest descent method to compute saddle points with one "de- 
scent direction" (such as a Mountain Pass solution). A proof of convergence of a variant 
of that algorithm was later given by Y. Li and J. Zhou in [10, 11]. To briefly describe it, 
let us fix a closed subspace E of Jtf? and <p a continuous £ ,± -peak selection, i.e. (p(u) is 
the location of a maximum of S on E®R + u := {e + ty \ e € E, t ^ 0} for any u E M'\E 
and <p is constant on £ © M. + u. As it will be convenient in the rest of the paper that <p is 
not solely defined on a unit sphere, we present here a slightly different version [21]. 

The authors are partially supported by a grant from the National Bank of Belgium and by the program "Qualita- 
tive study of solutions of variational elliptic partial differerential equations. Symmetries, bifurcations, singular- 
ities, multiplicity and numerics" of the FNRS, project 2.4.550. 10. F of the Fonds de la Recherche Fondamentale 
Collective. 
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Algorithm 1.1 (Mountain Pass Algorithm), (i) Choose uq e Ran 9, e > and n 

0; 

(ii) if\\ V(S'(m„)|| < e then stop; 
else compute 



u n+ \ = q>\u n -s„ 

for some s n €E S{u n ) C (0, +°°) where S(u n ) is a set of "admissible stepsizes " chosen 
so that at least the following inequality holds: 

£{u n+ {)~ S{u n ) < -±s n \\V£(u n )\\', 

(iii) let n <— n + 1 and go to step 2. 

Y. Li and J. Zhou proved that (u„) converges to a nontrivial critical point of $ 
up to a subsequence. The proof of convergence is performed in the space Jfif to ensure 
that the rate of convergence for the discretized problem does not deteriorate when the 
approximating subspace becomes finer. The original goal of the authors for introducing E 
was to try to obtain multiple critical points by taking E as the linear subspace generated by 
previously found solutions which the algorithm must try to avoid. The proof is performed 
in two steps. First, they show that s„ exists and that S decreases along (m„) h gn- This step 
relies on the following deformation lemma. 

Lemma 1.2. If (p is continuous, uq 6 Ran<p, uq ^ E, VS'(uq) 7^ 0, then there exists sq > 
such that 

Vs e M, 4<pK)) - *(«<>) < -y\\v*(uo)\\, 

where 

u s :=u -s . 

||v<r (mo) II 

The second step consists in proving, under some traditional assumptions on <p, that 
a subsequence of (u„) converges. For this, it is essential to show that the stepsize s n con- 
trols the distance between u„ and u n+ \ and that s„ is chosen in such a way that it is close 
to only when "mandated" by the functional. Let us remark that the choice of <p is very 
sensitive. Indeed, to seek sign-changing critical points, the modified mountain pass algo- 
rithm was introduced by J. M. Neuberger [14] (see also [7]). He considers algorithm 1.1 
above and only modifies the projection <p into a "sign-changing peak selection" (p^ which 
is a map defined from the set of sign-changing functions of Jif \ {0} to \ {0} such 
that, for any u, &(0pf(u)) > and (Pn(u) is a maximum of S on 0R + m~ where 

u + (x) := max{0, u(x)} and u~ (x) := min{0, u(x)}. Although in practice it appears to con- 
verge to a nontrivial sign-changing critical point, its convergence has yet to be formally 
proved. 

In this paper, (p(u) is allowed to pick up a maximum point of $ in an abstract cone 
C„ and we are interested in giving assumptions on C„ which imply the convergence of 
the mountain pass algorithm. This work is partly motivated by the article [17] where 
the authors define the notion of "natural constraints" to seek nontrivial critical points of 
functionals. Let us first make precise the peak selection <p that we use. We write intC 
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for the interior of C relative to spanC, the smaller closed subspace containing C, for the 
topology induced by Jf?. 

Definition 1.3. Let be an open subset of M '. We say that <p is a peak selection for 
(C u ) ues rf if (p is a map from srf to srf such that, for all u £ 

(i) C u is a closed cone pointed at 0; 

(ii) <p(u) £ intC„; 

(iii) for any v £ intC„, <p(v) = <p(m); 

(iv) (p(u) is a global maximum point of S on C„. 

Note that properties (ii) and (iii) imply that (p((p(u)) = <p{u). We write that d _L C„ 
if and only if d _L spanC„ for the inner product (-|). In section 2, we assume that <p is 
continuous and that (C„) verifies the following conditions: 

Vm £ .c/, « £ C u and 



3y £ (0, -), 38 £ (0, 1), Vmo £ Rancp, 3r > 0, Vm £ Ran(pnB(u ,r), 
Vd £ B(0,r), d J_ Q , C i7o+c/ nB( M0 , r) C Q + [1-5,1+ 5] A y {d) 
where A r (d) := {d' | = \\d\\ and Z{d',d) sC 7} and Z(d,d') := arccos(^^- 



(ACi) 
(AC 2 ) 

de- 



IMIKII 

notes the angle between two non-zero vectors d and d' (we set Ay(0) := {0}). This as- 
sumption, which speaks about the behavior of the cones under small deformations, is 
essential to prove a deformation Lemma in this generalized setting (see Lemma 2. 1). This 
lemma ensures the non-emptiness of the set S(uq) of admissible stepsizes at uq which we 
now define. For any uq £ Ran <p such that V#(«o) + 1 0, we set 



>o):={ 



= U>o 



u s := Uq — s- 



£ si and g(fp(u s )) -#0o) < — V^(«o)|| } 



|V<f(uo) 

for some value a > given by Lemma 2.1 and we require that s £ 5 (mo) := S*(uq) n 
[jSup5*(Mo),+oo). Other definitions of admissible stepsizes are possible provided they 
imply a local uniformity in the sense that s„ £ S(u n ) forces the stepsize s n not to be small 
when the gradient is not (see Lemma 2.4). The definition given above draws its inspiration 
from a paper [21] written by N. Tacheny and C. Troestler. 

To obtain the convergence of (u n ) up to a subsequence (see Theorem 2.14), we 
unfortunately need to replace (AC2) with the following stronger assumption: 
there exists a closed subspace E C J$? (possibly infinite dimensional) and 
a family of ^-vector fields : &f — > £ x , i = \,...,k, for some k £ N, 
such that for all u £ and i £ {1, . . . , fe}, 

(i) the family (|<(w)) =1 is orthonormal; 

(ii) Vv6V„, £/(a)[v]eV„, where V„ :=span{^i( M ),...,^( M )}; 

(iii) Vv £ intC„ n .t/, §(v) = $(«); 

(iv) («|^(")>^0; 

(v) 3r > 0, £;/ is bounded on {m | dist(«, Ran <p) < r}(~}£/. 

For m £ the cone C„ is defined as C M := £ © {£,■ f,^, (h) | ^ for all /}. 

Here, the notation dist(M, <9,e/) stands for inf { 1 1 u — v|| | v £ ds/}. Let us remark condi- 
tions (AC3) (i), (ii), (iv) and (v) are already present (albeit somehow implicitly for (iv)) 
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in [17] in the context of trivial ^-subbundles intead of cones. The additional condi- 
tion (iii) is equivalent to Vv G intC„ n srf, C v — C u . This condition is rather natural to 
require in view of property (iii) of the definition of peak selection. This "finite presenta- 
tion" of the cones is used in Lemma 2.8 to ensure that the stepsize s n controls the distance 
between u n +\ and u n . 

As a particular case of (AC3), let us mention that we can work with a family of 
continuous linear projectors (see Proposition 2.12). This case is an abstract formulation 
of the setting of [2, 3] where the convergence (up to a subsequence) of a mountain pass 
type algorithm for systems has been announced. 

In Section 2.3, we are interested in the convergence of the whole sequence generated 
by the Mountain Pass Algorithm. To that aim, we need to refine the definition of S* in 
order to control §{<$(u n — S ||v""i| )) for any < s < s„. 

In Section 3, we illustrate our method with two semi-linear problems. The first ap- 
plication takes its inspiration from a paper due to A. Szulkin and T. Weth [20] in which 
the authors study the following Schrodinger problem 

-Au(x) + V(x)u(x) = \u(x)\ p ~ 2 u(x), jc<e£2, 
u(x) = 0, x G <3£2, 

where V : £2 — > K is such that is in a spectral gap of -A + V and 2 < p < 2* := 
(+00 when N = 2). They are interested in the existence of non-zero solutions on an open 
bounded domain £2 C M. or on £2 = R N (in the latter case, V is assumed to be 1 -periodic 
in each jc,, i = l,...,N). Solutions to this equation are critical points of the indefinite 
functional 

S'-.J'f'^R-.u^- f (\Vu(x)\ 2 +V(x)u{x) 2 ) dx-- [ \u(x)\ p dx, (2) 
2 Jo p Jn 

where Jt? — Hq(Q.). The first proof of the existence of non-zero critical points for § 
when —A + V is not positive definite and £2 is an open bounded domain is due to P. H. 
Rabinowitz [18]. Recently, A. Szulkin and T. Weth proposed an alternative method [20] 
that also makes easier to deal with the lack of compactness that occurs when £2 = Mr. 
Denoting E the negative eigenspace of —A + V, they introduce the following nonlinear 
map 

<p : je\E -> Jf? : «H <p(m) 

where (p(u) is the point at which S reaches its maximum value on E ® E + m. They prove 
that minimizing § on Ran<p = \u € ,¥£\E | dS{u)[v\ = for v = u and any v € E\ 
yields a non-zero solution with least energy. Notice that, here, E is used to deal with the 
indefiniteness of the problem and not to compute multiple critical points as in the papers 
of J. Zhou & al. [10, 11]. We will prove that our algorithm converges for this problem. 
The numerical solutions that we obtain lead to some conjectures on the symmetries of 
ground state solutions. 
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The second application is based on a paper by B. Noris and G. Verzini [17]. The 
authors study the superlinear Schrodinger system 

{— Aui(x) = diF(u\ (x), ...,Uk(x)), x £ £2. 
, s „ 3 „ 1 = 1,...,*, (3) 

Ui(x) = 0, xe dQ., 

where k EN. They require that £2 C M. N is a bounded smooth domain and F € < ^ 2 (M.*; M). 
Note that the system — Auj = fijU? + UiY^j^i Pijtf where > and J3,j = jijj is a partic- 
ular case of (3). Such type of nonlinearities have been studied due to their applications to 
nonlinear optics and to Bose-Einstein condensation (see [6, 5, 8, 22]). Solutions to (3) are 
critical points of the functional 

: u = (u\ , • • • , ujA h ^ — / |Vi<(x)| 2 dx— f F(u)dx, (4) 
2 in Jn 

where Jf? = i/J (£2;R*). As already mentioned, B. Noris and G. Verzini [17] propose a 
general method of "natural constraints". Applied to the above problem, it goes as follows. 
Denote si :={«£ ffi \ U( ^ for every /}. To find a solution u — {u\, . . . , u^) of (3) 
with iij ^ for all i = 1, . . . ,k, they minimize £ on the constraint j¥ := | V/ — 

l,...,k, / Q |Vit,| 2 djc= f n diF(u)ujdx}. We will show that, under their assumptions, Ji 
Ran (p with <p being the peak selection 

(p '. si — y si \ u i — y argmax{i#(fiMi, . . . ,tk u k) \ h > 0,...,^ > 0}. 

Again, we prove that our algorithm converges for this problem and some numerical ex- 
periments are performed. 



2. Steepest descent method on varying cones 
2.1. Uniform deformation lemma 

Let S : ffi — > R be a ^ -functional defined on a Hilbert space and si an open subset 
of The following lemma is instrumental in proving the convergence of the algorithm. 

Lemma 2.1 (Uniform deformation lemma). Let <p : si — > si be a peak selection for 
(C u )ues/ an d uo € Ran<p foe smc/i that V£(uq) ^ 0. Assume that (p is continuous at uq and 
that (AC\ )-(ACt) hold. Then there exist sq > and ro > swcn f/zaf, /or any s € (0, So] an£ i! 
So G 5(w0) '"o) nRan<p, one has 

• V^(m )^0, 

• fij £ si where u s := fio — s ||v#(a°)|| anfl " 

• there exists some a > solely depending on J and 8 given in assumption (AC2) such 
that 

g (<p(fi t )) -£{u ) < -cw||V#(m )||- (5) 

Proof. Let «o £ Ran 9 C ^ and let us consider 7, 5 and r given by the assumption (AC2) 
for uq. Since ^ is open, there exists £1 > such that for any u E B(uq,Ei) and v G B(u,E\), 
one has «,v6^, V<f (a) 7^ 0, V#(v) ^ and «,v6 B(u ,r). 



6 



Ch. Grumiau, Ch. Troestler 



For any «eB(wo,£i), let d u := -V<? (k)/||V<?(w)||. Then 

VmgB(mo,£i), Vrf€A y (4), {VS(u)\d) sC -cos7||V(f( M )||. 

Let 7 := i cos 7 > 0. Taking £1 smaller if necessary, we may assume that 

V M ,vgB(m ,£i), VdeA r (4), (V<f(v)|d) < -f||V<?( M )||. 

Thus, on one hand, there exists £2 > such that, for any u G B(uo,£2), v G B(u,e2), 
d G A y (d u ) and <7 G (0, £2), 

(V#(v + arf)|d)<-7]|V<r( M )||. 

For any fio € B(uo,£i) HRan<p, v G Ca nfi(iio,£2), ^ G A r (c4 ) an d 17 < £ 2> tne mean 
value theorem implies there exists a a G (0, a) such that 

£{v + ad)-g{uo) ^£{v + ad)-g{v) (6) 
= (V£(v+6d)\od) 

<-fa\\V£(u )\\, (7) 

where the first inequality results from v G Ca and uo = (p(uo) is a global maximum of <§ 
onC l7o . 

On the other hand, by the continuity of <p at uq, there exist so G (0,r) and £3 G 
(0,min{r, ^£2}) such that, for any uo G B(uo,£i) an d s G [0,so]> one has <p(fio + *rfa ) G 
B(wo,min{r, 3 £2})- Let u s := uq +sda - If in addition uq G Ran<p, one has da _L spanQ,, 
(because So = (p(uo) G intCa, * s a l° ca l maximum) and therefore one deduces from as- 
sumption (AC2) that 

q>(u s ) G Cs s nS(M ,r) C C l7o + [1 - 8, 1 + S]A y (srf So ). 

Thus, q>(u s ) = v s +K s sd* for some v s G C [7o , G [1-5,1 + 5] and d* G A r (d fio ). So, 
possibly taking sq smaller, we get that K s s < ^£2 and v s = (p(u s ) — K s sd* G B(uo,£i)- 
Using equation (7), we conclude that 

<?(?(«,)) - #(fi ) = #(v s +^<) -<f (mo) -7(1 - S>||V<f (w )|| 

for any uq G B(mo,£3) nRancp and s G (0,so]- D 

Remark 2.2. • Equation (6) is the unique place we use that (p(u) is a global maximum 
of d? on C„. This assumption can be weakened by only requiring that the neighbor- 
hood on which (p(u) achieves the maximum of S is locally uniform w.r.t. u: 

Vmo G Ran<p, Bp > 0, Vm g Ran<pnfi(M ,p), £ (<?(")) = max «?(v). (8) 

v6C B nB(«,p) 

This assumption allows the existence of multiple maximums points in C u . It was not 
used in definition 1.3 for simplicity but also because, in the examples of section 3, 
(p(u) is a maximum on the whole C u . 
• Let us also note that, if we are just interested in the inequality (5) at uq (and not 
for all uo in a neighborhood of uo), we only need to require that (p(u) is a local 
maximum of <§ on C„. 
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• A careful reader may notice that we did not really use the fact that C„ is a cone 
pointed at 0. However, if (C u ) was just a family of sets satisfying {AC\), (AC2), (8) 
and the fact that <p(w) G intC„ in a locally uniform way: 

Vmo G Ran<p, Bp > 0, Vm G Ran<p nB(wo,p), B(u,p) n spanC„ C C„, (9) 

then the cone C„, defined as the closure of {t v \ t ^ and v G C M }, also satisfies (ACi ), 
(AC2), (8) and (p(u) G intC M . So very little is gained by not using cones, especially 
because they are the natural structures encountered in our examples. 

• As a consequence of the above deformation lemma, one can interpret Ran <p as 
somewhat a natural constraint for S in the sense of [17]. More precisely, it im- 
plies that if mo G Ran <p is a local minimum of S on Ran <p then uq is a critical point 
of £ on the whole space . 

2.2. Convergence up to a subsequence 

In this section, we first remark that it is possible to construct a sequence of stepsizes s„ 
such that the energy § decreases along the sequence (m„)„ g n generated by algorithm 1.1. 
In the following, without loss of generality, we can assume that V<S(u n ) 7^ for any n G N 
(otherwise the algorithm finds a critical point in a finite number of steps). 

Proposition 2.3. If s n > verifies inequality (5) given in Lemma 2.1 for any n G N, then 
the functional S decreases along the sequence (u n ) ne jq. 

Proof. As V ^= 0, s n is well-defined by Lemma 2.1. By construction, we have 

<g{u n+x ) - S{u n ) = S ((p(u n -s n ™ Y" „ ) ) - < -0M n ||V<?(v)|| < 0. 

So, <^(w„+i) < <£{u n ). □ 

As explained in the introduction, we now consider the sets S* (wo) and S(uq). The 
set S* (uq) is not empty as soon as uq is not a critical point of § (thanks to the deformation 
lemma). Concerning the set S(uq), it is not-empty once S is bounded from below on 
Ran 9, an assumption that we will later make (see Theorem 2.10). 

Lemma 2.4. If uq G Ran<p, VS (uq) and <p is continuous at uq, then there exists an 
open neighborhood V ofuQ and s* > such that S(u) C [s* , +00) for any u G V D Ran <p. 



Proof. By the uniform deformation Lemma 2.1, there exists sq > and ro > such that, 

WW 

||V#(h) 



for any < s ^ so and m G B(uQ,ro) n Ran <p, we have u s :~ u — s ,,Z £["|n G jz/, V<?(m) 7^ 



and 

<f (<P(m*)) -<?(m) < -CM||V<f(K)||. (10) 

In particular, for any m G B(uQ,ro) n Rancp, so G S*(u). It follows that C +°°). It 

suffices to take s* ^ so /2. □ 

Remark 2.5. To prove Lemma 2.4, let us remark that we could only use inequality (10) at 
u = uq for s = so fixed. Indeed, by continuity, it directly implies that so G S(u) for u close 
to uq. However, to obtain Lemma 2.4 for S(u) (see section 2.3) instead of S(u), the full 
strength of the deformation lemma is needed. 
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From now on, we have to require condition (ACj,). Let us first show it subsumes (AC2). 

Lemma 2.6. Let (|j)*_j be the family of vector fields given by (AC3) and assume (AC\) 
holds. Then 

k k 
Vm e<Vdl C u , £(« I &(n)) • %(u)[d\ =</-£(« I §/(«). 

i=l i=l 

Proof. For any m e the fact that »eC„C and that (|i)f = i ^ s an orthonormal basis 
of Vj, imply m = L( M I^!'( M ))^' ( M )- Differentiating in a direction <i £ Jf, yields 

</ = £<rf|&(«)> ^,(")+L("l^(")M>^(")+L("l^("))-4'(")M- 

i=l i=l i=l 

If of _L spanC,,, the first term vanishes. This completes the proof. □ 
Proposition 2.7. Properties (AC\) and (AC3) imply (AC2). 

Proof. Let 8 £ (0, 1) (property (AC2) will be satisfied whatever value is chosen). Simple 
geometrical considerations show that there exists a 7 £ (0, 7t/2) such that 

B(d,8\\d\\) C [l-5,l+5]A y (</). 

Let «o € Ran<p. As mo £ intC„ , there exist a > such that (wo||('(mo)) > 0£ for all z. 
Using the continuity of §j and at mo, we can choose r sufficiently small and a M > 
(depending only on mo) so that, for all u £ B(uo,r) and all i = l,...,k, 

(a|&(a))>a, ||6(h) -&(fio)||<e, |||/(k)KM, and ||#( H ) - |/(m )|| < £, 

where e > is a constant depending only on 8 and uq (to be chosen later). 

Let mo € B(uQ,r) and c/ g B(0, r) such that d _L C« . Let w £ Cu +d P\B(uo,r). One 
can write w = e + ££f|f(z<o + d) for some e £ E and f,- ^ 0. Let us start by noticing that 
tj = (w\%i(uo + d)} . Therefore, recalling that ||^,|| = 1, one deduces 

(fio||/(fio))| sc \(w-uo I §(ua+d)) \ + \(uo\ %i(u +d) -§(wo))| 

5$ 1 1 w — w 1 1 + ||moIIII4/(«o + ^) -|i("o)|| 

<2r+(||no|| + r)e. (11) 

Provided that e and r are chosen small enough, one can assume that 2r+ (\\uo\\ +r)e ^ 
a/3. In particular, this implies f, > 2a/3 > 0. 

Using the integral form of the mean value theorem, we get 

k k pi k 

w = e + £ ^(wo + d) = e + £ r 4'(«o) + / £ f i 1/ ( fi o + *0 M ds - ( 12 > 

i=i 1=1 Jo ;=i 

The third term can be rewritten as follows: 

k k 

£ («o | |;(«o) ) 1/ (mo) M + L ( f '' ~ ("0 ||;("o) )) 1/ («o) [d] 
t=i t=i 

,1 * 



/ £ *i (1/ ("0 + «/) [d] - 1/ (So) [rf] ) ds =: rfi + c/ 2 + * • 

J ;_i 
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Using Lemma 2.6 on d\, one can write equation (12) as 

k 

w = e + Y,(ti-{uo\%(uo)[d}))Zi(uo)+d+d2+d 3 . 
i=i 

Since | (So|^/(fio)[c/]) | ^ ||£o||||€/0<o)|| \\d\\ ^ (|| M o|| + r)Mr, we can assume r was chosen 
small enough so that this is smaller that a/3. Recalling that f, > 2a/ 3, one sees that the 
coefficients of ^,(fio) are positive and therefore e + L(f/— (mo|^/(5o)[^])) <^'("o) & Ca Q . 

The proof is complete if we show that d + d2 + di G B(d,8\\d\\). Using (11), we 
deduce |f,-| ^ \ti — (uo\%i(uo)) \ + ||«o|| ^ 2r+ (||i<o|| +r)(e + 1). Thus the following esti- 
mates 

k 

( fi o||i(«o))| U!(uo)\\\\d\\^ *(2r+(||« || +r)e)M\\d\l 

i=i 

k 

||<fc||<J>| sup ||4/(Mo+^)-^("o)||||« f Kfc(2r+(|| Mo ||+r)(e + l))£|| fl f||, 

;=i se[o,i] 

show that ^ jSII^H, i = 2,3, provided that the constants e and r were chosen small 
enough. □ 

Lemma 2.8 is the second key element to prove the convergence up to a subsequence. 

Lemma 2.8. Let <p be a peak selection for (C u ) ues / which verifies conditions (AC\) and 
(AC3). Let (wn)neN an d (s n )„ e jq be given by the generalized MPA (algorithm 1.1) with 
s„ G S(u„) for all n. Let us assume that <p is continuous, Ran<p C si ', and either 

3ti,...,^G (0, +00), dist({ Y,z&(u) KGRanp},^) > 0, (13a) 



1=1 

V(v„) C Ran<p, (<o(y n ) \ is bounded from above => (v„) is bounded 
and dimE < °°. 



(13b) 



IfY.t=o s n < +°° tnen ( M «)neN converges in si '. 

Proof. Let k be given by the assumption (AC3). For i = l,...,k, set v,„ := and 
d„ := — ||v#(''"]|| ■ Let r be given by assumption (AC3) (v) and A", be a bound for Denote 
K := max, = i ^Ki. 

By assumption (AC3) and as <p(M„ +s„d n ) G intC„ n+4ni /„, we have 

Vi,n+1 = 4;(<p(m„ +.s„rf„)) = t,i{u n +s n d n ) 
for any n G N. Let n* be large enough so that s„ < r. Thus, for all n ^ «*, 

||Vi,n+l-Vi )f! || ^^||5 n rf„|| = #,S„. (14) 

Since £s„ < +°°, it follows that for any i = l,...,k, (v,„)„ g n is a Cauchy sequence and 
therefore converges to, say, v !)0 o. 

Let us assume (13a) holds. Consider v„ :— YX=\ ti v i,n = L?=i Ti^i(u„). It converges 
and its limits belongs to sf. Since (p(v„) — (p(u n ) = u n and <p is continuous, the sequence 
(«n)neN converges. Its limit lies in Ran <p and thus in si ', 
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If on the other hand (13b) holds, the fact that the sequence (<£*(«,,)) is decreasing 
implies that (u n ) is bounded. Let (u' n ) be a subsequence of (u„). Since u' n € C u t , one can 
write u' n = e' n +Ef=i f ,><^( M «) f° r some e' n £ E and f ( ' n e (0,+°°). As (t^) is bounded, so 
are (ej,) and |f/J = |HI&( M «))I < IKII- So ' U P t0 subsequences, (e' n )„ m and (f,'„)„ GN 
converge to, say, e^, and t' ioo . Thus — > u'^ := e'^ + £fj qD v,-,»>. Thanks to Ran<p C 
mL € .e/. But then the continuity of and 9 imply 

v<,» = |,-(0 and = q>(u' n ) <p(iQ- (15) 

If the same reasoning is performed with another subsequence («"), (15) implies that 
^i( M 'e») — <=i'( M ~) an d therefore, in view of definition 1,3 (iii), (p(uL) — As the 

limit does not depend on the subsequence, the whole sequence (u n ) converges in sd '. □ 

Remark 2.9. If we wanted to seek sign-changing solutions using the cones C u := R + u + © 
M + m~ (as explained in the Introduction), then we would not be able to remove the projec- 
tion factors in the above computation of v,.„ + 1 . This sheds some light on the difficulty of 
proving the convergence of the Modified Mountain Pass Algorithm [14]. 

Theorem 2.10. Assume <p : srf — > is a continuous peak selection s.t. Ran<p C &/ 
and the cones (C u ) ueji / verify the conditions (AC\), (AC3) and (13a) or (13b). Suppose 
moreover that <§ 6 < j/s 1 (J%?;13L) satisfies the Palais-Smale condition in Ran<p and that 
inf„ e R an( p §(u) > —00. Then the sequence (u n ) ne f$ given by the generalized mountain pass 
algorithm 1.1 possesses a subsequence converging to a critical point of § in Ran (p. In 
addition, all limit points of (M„)neN are critical points of S. 

Proof. Let us start by showing that (V<?(« n )) „ converges to zero up to a subsequence. 
If not, we could assume there exist 8 > and «o s N such that, for any n ^ «o, || V#(w„)|| > 
8. Then, for any n ^ «o, the deformation lemma 2. 1 implies 

&(un+i) —& (u n ) sS -as n 8. 

Thus, summing up, 

+00 +00 
\xm^g(u n )-<g(u no )= s "- 

As the left-hand side is a real number (£" is bounded from below on Rancp and de- 
creasing along (m„)„ g n), we have Y^=o s " < +°°- So, by Lemma 2.8, u„ — > u* e &/ and 
||V<f (w*)|| > 8. By continuity of <p at u* G ^f, we obtain <p(m*) = «* and, so, u* e Ran<p. 
By Lemma 2.4, there exists a neighborhood V of u* and s* > such that 5(m) C [5* , +°°) 
for any u € V. Consequently, there exists «o such that, for any n ^ no, s n 5* whence 
^n=o s " does not conver g e ^ which is a contradiction. 

In conclusion, there exists a subsequence (u n ,)k e f$ of (w n )neN s.t. ||V(?(K n .)|| — > 
when — > +00. As <» satisfies the Palais-Smale condition, (w nj ,)fe=N possesses a subse- 
quence converging to a critical point of S . 

Concerning the second statement of the theorem, the argument is very similar. Let 
( u n k )keN be a convergent subsequence and assume on the contrary that u := lim^,*, w n . 
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is not a critical point of £ . In that case, on one hand, there exists 8 > and k\ £ N such 
that, for any k ^ k\, \\V<S(u nk )\\ ^ 8. By Lemma 2.1, we have 

Vk^k u S{u nk+l )-<g{u nk ) ^~a8s„ k . 

On the other hand, as u £ Rancp, we have by Lemma 2.4 that 

3 s* > 0, 3k 2 £N,yk^k 2 , s„ ^ s* . 

So, for large k, §(u nk+{ ) — §{u nk ) ^ — "8s*, which is a contradiction because n6N 
is a convergent sequence. □ 

Remark 2.1 1. By previous remarks 2.2 and 2.5, we conclude that we could get the con- 
vergence up to a subsequence using the equation (5) only at mo- Thus, the uniform form of 
Lemma 2.1 is not required (and we could consider that (p(u) is a local maximum of <p on 
C„ instead of a global maximum). Nevertheless, we have kept the uniform setting along 
the paper as it will be required in Section 2.3. 

The following special case of (AC3) is important for the applications. 

There exist a closed subspace E C Jff (possibly infinite dimensional) and 
linear continuous projectors P, : Jf —> E- 1 , i = l,...,k, for some k £ N, 
such that 

t Vh€ J4? , Pi(u) _L Pj(u) whenever i 7^ jf; 
• E®l$ =1 RaaPi = Jf. 
For all u £ Jf , set C„ := E © {E^(w) | f; > for all 2} . 



(AC 4 ) 



Let us now sketch the proof that (AC4) implies both {AC\) and (AC3). Consider 



s/:={u£Jtf\ Pi{u) ^ for all i} and |,-(k) := 4^. Clearly §1, . . . , % k are ^ functions 



on s/. Moreover e + £f,Pi(«) €intC„ifand only if all f,- > 0. Since w = Pe(u) +Lf=i 
where denotes the orthogonal projection on E, one has u £ intC M . Given the definitions 
of s/ and points (i), (iii) and (iv) of (AC3) are straightforward. A simple computa- 
tion shows that §/(w) [££/P/(h)1 is a multiple of P,-(m) whence (ii) follows. Finally, as 
||!/(k)|| = 0{\/\\Pi{u)\\), (v) will hold provided ||P,(m)|| is bounded away from when 
u £ Ran (p. Note that this latter condition also ensures that Rancp C si ', Remark that 
these cones satisfy property (13a) with Tj = • • • = = 1 because dist(£ §/( m )j — 
rmnj\\Pj(Ui(u))\\ = 1. 

Thus, as a corollary of Theorem 2. 10, we get the following proposition. It can be 
thought as an abstract version of the convergence results in [3, 10, 1 1]. 

Proposition 2.12. Let us consider <p : si — > si a continuous peak selection with the 
cones (C„)„ e ^/ being given by condition (AC4) and sf := {u £ Jj? \ Pi(u) 7^ Ofor all i = 
1, . . . ,k}. Assume that inf l(G R an( p||P(M)|J > Ofor all i — l,...,k, that S £ ^ (Jf ;R) sat- 
isfies the Palais-Smale condition in Ran <p and that inf„ e R an( p <S(u) > — °°. Then the se- 
quence (w K )„ e N given by the generalized mountain pass algorithm 1.1 possesses a subse- 
quence converging to a critical point of S in Ran <p. In addition, all limit points o/(w M )„ e pj 
are critical points of S. 
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In Theorem 2.10, the Palais-Smale condition is required. For the particular case of 
= H (Mr), this condition does not generally hold as mass may be lost at infinity. 
Fortunately, H 1 (M ) respects the following compactness condition (see for example the 
paper [12] for a proof): for any bounded sequence (m„)„ g n Q H l (M n ) staying away from 
zero, there exists (x n ) ne ^ C Z N such that (u n (- +x n y\ weakly converges up to a subse- 
quence to a non-zero function. This is enough to get the convergence up to a subsequence. 

Proposition 2.13. Assume the hypotheses of Theorem 2.10 hold, except for the Palais- 
Smale condition. Let ,34? := H (K ) and (u n ) ne ^ be the sequence given by the Moun- 
tain Pass Algorithm 1.1. If for any u € ,3i? and x € Z , S(u(- + x)j = §(u) and if 
ffl — > Jf? : u i — y VS(u) is continuous for the weak topology on -3^, then there exists a 
sequence (x„)„ G pj C Z N such that («„(• +*n)) nGN weakly converges up to a subsequence 
to a nontrivial critical point of §. 

Proof. We will only briefly sketch the proof. As (m„)„ g n is bounded in H 1 (R N ) and stays 
away from 0, the compactness condition recalled above implies that there exists a se- 
quence (x„) mG n C Z n such that u(-+x n ) weakly converges, up to a subsequence, to u* ^ 0. 
Intuitively, the translations "bring back" some mass that u„ may loose at infinity. 

Using the translation invariance of S, the corresponding equivariance of V<f and the 
weak continuity of VS 1 , we conclude that u* is a critical point of S . □ 

2.3. Convergence of the whole sequence 

In this section, we refine the stepsize used previously to get the convergence of the whole 
sequence generated by algorithm 1.1. We require that the stepsize s n £ S(uq) : — S*(uq) D 
Q sup 5* (mq), +°°) where 




Using the deformation lemma 2.1, we get that S(u„) ^ as long as u„ is not a critical 
point. Moreover, working as previously, we get results 2.4, 2.8 and 2.10 for this new 
choice of stepsizes. Let us remark that, this time, we really need that inequality (5) is valid 
in a neighborhood of mo to get Lemma 2.4. This new stepsize will allow us to control 
the energy for any < s ^ sq. Under a "localization" assumption, we now prove that 
the whole sequence (m„)„ g j*j given by the mountain pass algorithm 1.1 converges to a 
nontrivial critical point of S '. 

Theorem 2.14. Assume that u is the unique critical point of S in the ball B(u, 8) for some 
8 > 0. Under the same assumptions as those of Theorem 2.10, if there exists n* € N such 
that <§(uy?) < a := in{ v edB(u,8)r\Ran<p <*( v ) an d u n * € B(u,8) then the sequence (m„)„ g n 
produced by algorithm 1.1 with stepsizes s n € S(u n ) converges to u. 

Proof. For any m > n* , we claim that u m G B(u, 8). If not, as u n * 6 B(u, 8), there exists 
m ^ n* such that u,„ G B(u,8) and u m+ \ — (p(u m — s m pfg^" 1 ^ ) ^ B(u,8), with s m e 
S(u m ). By continuity, there exists < s ^ s m such that <p (u m — s ||vff(""'j|i ) € dB(u, 8)0 
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Ran (p. This is a contradiction because, by the definition of s m and as § is decreasing along 
(m„)„ g n, we have a < S (<p(« m - s ||v^")|| )) ^ ^(«m) ^ &(u„*) < a. 

As // is the unique critical point in B(u,8), by Theorem 2.10, u is the unique accu- 
mulation point of (w,,)neN- So, m„ converges to u. □ 



3. Applications 

3.1. Application to Indefinite Problems 

For problem (1), the energy functional S given by (2) is defined on 3? := Hq (Q.). Let us 
denote the decomposition Jt? = Jtf^ Jf corresponding to the spectral decomposi- 
tion of — A + V with respect to the positive and negative part of the spectrum. For any u, we 
let u^> € 3i?{> and £ Ji?( + > be the unique elements such that u — u^"> + u( + K Let 
us remark that the case J4?(~> — {0} corresponds the traditional mountain pass algorithm 
with a positive definite linear operator. 

We choose the following peak selection. Let srf := and, for any u G sd ', 

let C„ be the cone C u := H 0l + M = (_) © The peak selection (j? for (C„) 
is the map 

q> : srf — > si : u i-> <p(w) 
such that, for all m g jz/, ^)(m) maximizes S on C„. To prove that <p is continuous, we 
refer to the original paper [20]. Is is easy to check that these cones verify (AC\). Indeed it 
suffices to consider E = Jff^~\ k=l and Pi : J$? — > £ x , the orthogonal projection on £ x . 

To apply Proposition 2.12, we need to verify the following assumptions on $: on a 
bounded domain i2, 

(i) it is standard to show that # G < i x ; 

(ii) S verifies the Palais-Smale condition on Rancp (see [20]); 

(iii) inf„ e R an( p S(u) > actually $ is bounded from below by on Rancp, see [20]; 

(iv) does not belong to Ran Pi o <p: it comes from the fact that is a strict local mini- 
mum of § on £ x = JfW (see [20]). 

In conclusion, Proposition 2. 12 applies and gives the convergence up to a subsequence of 
the sequence (u n ) generated by generalized mountain pass algorithm 1.1 for this indefinite 
problem provided that the domain Q, is bounded. 

Let us now sketch what happens about the convergence up to a subsequence when 
Q. = ~R N . As (<?(wn))nGN is decreasing (see 2.3) and is bounded away from zero, we have 
that (m„)„ g n is bounded and stays away from zero in H (Mr) (see [20]). On the other 
hand, V is assumed to be 1 -periodic, thus $(u(- +xj) =£{u) for any u € Jtf' andx £ I* N . 
It is not difficult to check that V<f is weakly continuous. Thus, Theorem 2.13 asserts 
that, if (u n ) is the sequence generated by the MPA, there exists a sequence of translations 
(x n ) C 1 N such that («„(• +x„)) nG pj weakly converges, up to a subsequence, to a nontrivial 
critical point u* of S . Moreover, if S(u n ) — > inf M eRan<p <S(u), then it can be proved that the 
above convergence is strong. The idea is that, if it does not converge strongly, some mass 
is lost at infinity. At the limit, this mass will take away a quantity of energy greater or 
equal to inf„ eRan( p §{u) > 0, a contradiction. 
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Numerical experiments. Let us start by giving some details on the computation of vari- 
ous objects intervening in the MPA. Functions in ffi will be approximated using P l -finite 
elements on a Delaunay triangulation of Q. generated by Triangle [19]. The matrix of the 
quadratic form (u\ , U2) f a Viq V«2 is readily evaluated on the finite elements basis. For 
(u i , U2 ) i-> V (x) u i «2 dx and the various integrals involving u to a power, a quadratic in- 
tegration formula on each triangle is used. The gradient g := V<?(v) is computed in the 
usual way: the function g g ffl is the solution of the linear system of equations V(f> g ffl , 
(g\(p),yf — d<f(v)[<p]. In practice, the peak selection <p must be evaluated with great accu- 
racy to obtain satisfying results. For this, we use a limited-memory quasi-Newton code for 
bound-constrained optimization [13]. The program stops when the gradient of the energy 
functional at the approximation has a norm less than 10~ 4 . 

As an illustration, we consider Q. = (0, 1 ) 2 , V £ K constant and p = 4. Let us remark 
that is then formed by eigenfunctions of —A + V with negative eigenvalues. In 

dimension 2, the eigenvalues A, of —A on the square (0,1) 2 with zero Dirichlet boundary 
conditions are given by n 2 (n 2 + m 2 ) with n,m = 1,2, . . . The related eigenfunctions are 
given by sin(n7Dc) sin(mny). We get Ai = 2n 2 w 19.76, A2 = A3 = 5n 2 « 49.48 (a double 
eigenvalue), A 4 = $% 2 w 78.95, A 5 = A 6 = 10^ 2 w 98.69,... 

Figure 1 depicts four non-zero solutions approximated by the algorithm 1.1 for four 
different values of V. The algorithm was always started from uo(x,y) := xy(x~ l)(y — 1). 
The graphs on the left-hand side are given for the values V = (dimJ^^ ) = 0) and 
— A2 < V = —21 < —A] (dimJ^~) = 1). The graphs on the right-hand side are given for 
-A4 < V = -50 < -A 3 (dim Jf? (-) = 3) and -A 5 < V = -80 < -A 4 (dim^'-' = 4). 
In Table 1, we present some characteristics of the solutions. 



V 


||V#|| 




# of steps 


£{u) 





6.0-10" 


-5 


7 


37.89 


-21 


6.4- 10- 


-5 


48 


70.43 


-50 


5.3-10" 


-5 


113 


91.42 


-80 


6.5 ■ 10" 


-5 


44 


35.06 



TABLE 1 . Characteristics of approximate solutions to an indefinite problem. 

For V = 0, we remark that the approximation is even w.r.t. any symmetry of the 
square and is positive. It was expected and it is actually already known in this case (i.e. 
for the problem —Am = \u\ p ~ 2 u) that ground state solutions have the same symmetries as 
the first eigenfunctions of —A (see [9, 1]). 

For V = —21, the approximation has two nodal domains and a diagonal as nodal 
line. It seems to respect the symmetries of a second eigenfunction of —A. It can be ex- 
plained as follows. When V — 0, it is proved [1] that, for p close to 2, least energy nodal 
solutions have the same symmetries as their projections on the second eigenspace of —A. 
On the square, it is even conjectured that the projection must be a function odd w.r.t. 
a diagonal. In view of the bifurcation diagrams computed by J. M. Neuberger [15, 16], 
the least energy nodal solution for V G (— Ai,0] becomes the solution with lowest energy 
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y=0 V = -50 




Figure 1 . MPA solutions for an indefinite problem on a square 

when V E (— A2, — Ai] and no bifurcation happens along the way. So it is reasonable (and 
this is supported by the bifurcation diagrams) that they keep the same symmetries along 
the whole branch. 

We also observe that, for V = —50 (resp. —80), the approximation seems to respect 
the symmetries of (and has the "same form" as) a fourth (resp. fifth) eigenfunction of —A. 
Their number of bumps corresponds to their Morse index (dim Jif (~> + 1). 

All those considerations support the conjecture that if — A„ < V < — A„_i then, at 
least for p small enough, ground state solutions respect the symmetries of a n th eigenfunc- 
tion of —A. 

3.2. Application to Systems 

In this section we will perform numerical experiments for the system (3). The corre- 
sponding energy functional (4) is defined on — ffi(i2,R fc ) endowed with the norm 
II "II 2 = in I ^ M 1 2 = Hila\^ u i\ 2 ^ x - 1° [17]> B- Noris and G. Verzini prove that the min- 
imization of § on Jif := \u € srf | V/ = l,...,k, J n \Vui\ 2 dx = f a diF(u)ui6xj, where 
£/ :={u£ Jrf? | K; ^ for every /}, yields a solution u = (u\, . . . = Y. u i e r with w; ^ 
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for all i = 1 , . . . , k provided that the following assumptions are satisfied: there exist p £ 
(2,2*), O > and 8 > such that, for any u, X £ R k , one has 

(0 LiMjF( U )\ < C f \u\p-\ E,|^F(«)| < C f \u\p- 1 and |F(«)| < C F |«| P ; 

(ii) LijdfjF^XiUiXjUj - (l+8)Z i d i F(u)Xlu i > 0; 

(iii) for every i there exists S, > such that <9,F(S,e,) > 0; 

(iv) diF(u)iii ^ diF(ujei)ui for every z. 

The first three assumptions are traditional in the framework of variational methods. The 
last one insures m, =^ for all i. In this section, we will use the Mountain Pass Algo- 
rithm 1 . 1 with the following peak selection. For any u = (hi, . . . , «&) G si , we consider the 
cone C u :={(t\U\,..., t^u^) \ ti ^ for all i = 1, . . . ,k}. The peak selection <p for (C u ) uea / 
is the map 

q> : si — > si : u >-» <p(w) 

such that <p(m) maximizes on C„. Under the additional hypothesis that Y,i^iF(u)ui ^ 0, 
the second assumption plays the role of the traditional super-quadraticity and implies that 
<p is well-defined as apeak selection. In fact, if u g si verifies &§(u)\(\\u\, . . . ,A&i<&)] = 

for any (Ai , , A^) € M. k then u is a strict local maximum of S on C u . It implies the 

uniqueness of the global maximum of S on C„. Moreover, <p is continuous. 

To see that assumption (AC4) is satisfied, it suffices to take E = {0} and, for i = 
l,...,k,Pi(u) =Pi((ui, . . . := Ujej i.e., P; is the projection on the / th component of u. 
Finally, let us quickly run through the assumptions of Proposition 2.12: 

(i) it is standard to show that S € ^f 1 ; 

(ii) S verifies the Palais-Smale condition on Rancp (see [17]); 

(iii) inf„ e R an( |, S(u) > actually £ is bounded from below by on Rancp (see [17]); 

(iv) dist(Ran<p,<?.K/) > (see [17]); 

In conclusion, Proposition 2.12 applies and gives the convergence, up to a subsequence, 
of the sequence (m„) generated by the Mountain Pass Algorithm 1.1. 

Numerical experiments. For the numerical experiments, we will consider the following 
particular case of equation (3): 

-Aui{x) = Hiii] + u t £ PijUj, x e £2, 

m i=l,...,k, (16) 

ui(x) = 0, x e d £2, 

where /},• = j3; j and £2 is a bounded domain of R 2 . This system is modeling a competition 
between A: populations. We will focus on the case Q. = (0, l) 2 and k = 2. In this setting, 
the assumptions (i)-(iv) stated above boild down to 

jlti>0, H2>0, and - y/jMVa < J3i,2 < 0. (17) 

Let us remark that the condition Y,i^iF(u)uj ^ discussed in the previous section is also 
verified in this range. 

Let us now give the outcome of the algorithm for various choices of (/ii,^a,j3i 2). 
The MPA will always start with the function uq = («o.1j m o,2) G where uoi(x,y) = 
u 02(x,y) — xy(l — x)(l —y) and stops when the norm of the gradient is less than 10~ 4 . 
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First we choose (/Ji,/i2>j3i,2) = (1,4, — 1). The numerical solution (mi, 1*2) is de- 
picted on Figure 2 and some characteristics are given in Table 2. In this case, the assump- 
tions (17) are satisfied so the fact that the algorithm converges to a solution (ui,ui) with 
u\ > and > is expected. Notice also that the solutions u\ and U2 are even w.r.t. axes 
of symmetry of the square. 

As second choice, we consider ,fl2,Pi,i) = (1:4,0.5). The MPA solution (111,112) 
is depicted on Figure 3 and some characteristics are given in the second row of Table 2. 
Despite the fact that the assumptions (17) are not satisfied anymore, the solution is sim- 
ilar to the found in the first case. If we enlarge j3i 2 further and choose ,/^2, j3i.2) — 
(1,4, 1.2), the algorithm still converges (see the third row of Table 2) but this time, the 
second component vanishes (see Figure 4). What happens is that, at the very first step, 
112 = and then the MPA essentially proceeds as if the system was only consisting in the 
first equation. 




Figure 2. MPA solution for the system with ^2,^1,2) = (1,4,-1). 




Figure 3. MPA solution for the system with (1x1,1x2,^1,2) = (1,4,0.5). 
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Figure 4. MPA solution for the system with (/ii, ^2,^1,2) = (1,4, 1.2). 



(Ml, M2, £1,2) 


||V*(«)|| 


# steps 


g(u) 


max mi 


maxM2 


(1,4,-1) 


7.9 -lO" 5 


11 


88.4 


8.6 


5.4 


(1,4,0.5) 


5.4 -lO" 5 


11 


40.4 


6.4 


2.4 



(1,4,1.2) 5.2 -10- 5 11 39.9 6.6 0.0 
Table 2. Characteristics of the solution to system (16) 
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